SPECIFICATION 
ADAPTIVE FILTER TO REDUCE MULTIPATH 

Cross Reference to Related Applications 

This method of reducing multipath avoids the cross term noise produced by other 
methods [Refs. 1 through 10]. These other methods correlate (cross correlate or auto 
correlate) two versions of the signal each containing the direct path and multipath 
reflections. A correlation peak between the direct path and a reflection gives a 
measurement of the relative delay between the direct path and the reflected path and the 
amplitude and phase of the reflection. As the relative delay between the two signals is 
moved from this peak, other correlation peaks are produced when there is correlation 
between the reflections, i.e., when the relative delay equals the difference of the delays 
for two objects. This produces noise which interferes with the measurement since there 
are many possible combinations. When one delays one of the signals and subtracts it from 
the other with the correct amplitude and phase, he tends to cancel the multipath. The 
other reflections, however, are dragged along with it and produce noise. These noises are 
avoided by the invention described in this disclosure. 

Statement Regarding Federally Sponsored R&D 
No Federal funds were used to conceive or develop this invention. 

Reference to Appendix Containing Computer Program 

The Appendix contains listings of computer programs written in the APL language and 
a sample computer run. The programs embody the basic features of this invention. The 
computer run is a demonstration of its performance on simulated data. 

The program NORM generates simulated data. The program DELRES calculates the 
delay residuals. The program NLS minimizes the residuals defined in DELRES by 
nonlinear least squares. 

The program NLS is first edited so as to output the residuals as RR. The program 
NORM is used to generate 100000 mean zero variance 1 random numbers in X5. This 
represents a 50000 Hz broadband signal sampled at 100000 samples per second with 
r.m.s. (root mean square) of 1. Two multipaths are simulated both of 0.1 magnitude. One 
is delayed 1 sample; the other is delayed 3 samples. The result is put into XX5. The 
signal to multipath noise ratio is calculated to be 17 dB. The NLS program is applied to 
XX5 and the two multipath noises are canceled to a level of 5 1 dB signal to multipath 
noise ratio. 
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Background 

The direct path of a radio signal from transmitter to receiver is frequently interfered 
with by reflections of the signal from stationary and moving objects. This is called 
multipath noise. This invention utilizes a new adaptive filter technique to reduce 
multipath noise. This method of reducing multipath avoids the cross term noise produced 
by other methods [Refs. 1 through 10]. These other methods correlate (cross correlate or 
auto correlate) two versions of the signal each containing the direct path and multipath 
reflections. A correlation peak between the direct path and a reflection gives a 
measurement of the relative delay between the direct path and the reflected path and the 
amplitude and phase of the reflection. As the relative delay between the two signals is 
moved from this peak, other correlation peaks are produced when there is correlation 
between the reflections. This produces noise which interferes with the measurement. 
When one delays one of the signals and subtracts it from the other with the correct 
amplitude and phase, he tends to cancel the multipath. The other reflections, however, are 
dragged along with it and produce noise. These noises are avoided by the invention 
described in this disclosure. 

Brief Summary 

A nonlinear least squares method measures the delay, Doppler shift and amplitude of 
the multipath due to each object and subtracts a very accurate reconstruction of each 
multipath signal from the noisy signal. If an object is a target, its range, range rate and 
magnitude is got from the measured multipath delay, Doppler shift and amplitude. 
Position and velocity of the target can be obtained by geometric triangulation with 
multiple transmitters. Target angle can be measured by the relative phase of the 
corresponding filter coefficients across multiple antennas. The system can be used on a 
surveillance aircraft to cancel ground reflections and measure targets. 

Brief Description of Drawing 

The Figure is a block diagram illustrating the adaptive filter of the invention to reduce 
multipath in a radio transmitted signal wherein the multipath is caused by reflections 
from fixed and moving objects, interferers, or targets. 

The method works best on white signals; therefore, prewhitening and post unwhitening 
filters are employed. 

The differential delay and Doppler shift for each target is measured by the filter weight 
x mn which corresponds to the differential delay nDand Doppler shift mf caused by that 
target. The magnitude of the filter weight is a measure of the strength of the target. 
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The quantity labeled "Residual" in the Figure is the received signal minus the result of 
application of the filter weights w, , w 2 , w 3 , x 1 , , x n , x l3 • • • , x 2] , x 22 , x 23 • • • jc 31 , x 32 , x 33 ■ - • .A 
nonlinear least squares method is used to pick the filter weights which minimize the mean 
squares of the "Residual". The result is a clean version of the original signal which is the 
received signal minus the multipath reflections. 

The filter weight x mn corresponding to a moving target is a measure of the differential 
range nDc and range rate mfc where c is the speed of light. The system, with a common 
antenna, can be duplicated to receive two signals each from each of two widely separated 
transmitters. Geometric triangulation can be used to measure the two dimensional 
position and velocity of the target. Triplication can be used for three dimensional position 
and velocity. 

Another use of the system is to detect the angle of moving targets. If multiple antennas 
are provided, each connected with a system like that shown in the Figure, the angle of a 
moving target causing a multipath reflection to the receiving antennas can be measured 
by using the relative phase of the corresponding delay and Doppler complex coefficients 
across the several antennas . 

The receiving antennas, receivers and processing system can be placed in a surveillance 
aircraft the position and velocity of which is obtained by an accurate navigation system 
such as GPS. Objects on the ground, interferers and targets causing multipath reflections 
can be processed by the system. Target position and velocity can be obtained by adding 
the position and velocity of the surveillance aircraft to the measured position and velocity 
of the target. 



Detailed Description 

This method of reducing multipath avoids the cross term noise produced by other 
methods [Refs. 1 through 10]. These other methods correlate (cross correlate or auto 
correlate) two versions of the signal each containing the direct path and multipath 
reflections. A correlation peak between the direct path and a reflection gives a 
measurement of the relative delay between the direct path and the reflected path and the 
amplitude and phase of the reflection. As the relative delay between the two signals is 
moved from this peak, other correlation peaks are produced when there is correlation 
between the reflections. This produces noise which interferes with the measurement. 
When one delays one of the signals and subtracts it from the other with the correct 
amplitude and phase, he tends to cancel the multipath. The other reflections, however, are 
dragged along with it and produce noise. These noises are avoided by the invention 
described in this disclosure. 
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The Figure is a block diagram illustrating the adaptive filter of the invention to reduce 
multipath in a radio transmitted signal wherein the multipath is caused by reflections 
from fixed and moving objects, interferers, or targets. 

The method works best on white signals; therefore, prewhitening and post unwhitening 
filters are employed. 

The differential delay and Doppler shift for each target is measured by the filter weight 
x mn which corresponds to the differential delay nDand Doppler shift mf caused by that 
target. The magnitude of the filter weight is a measure of the strength of the target. 

The quantity labeled "Residual" in the Figure is the received signal minus the result of 
application of the filter weights w, , w 2 , w 3 , x x , , x 12 , x 13 • • • , x 21 , x 22 , x 23 - • x 31 , x 32 , x 33 • • • .A 
nonlinear least squares method is used to pick the filter weights which minimize the mean 
squares of the "Residual". The result is a clean version of the original signal which is the 
received signal minus the multipath reflections. 

Refer to the Figure. The signal is modulated and amplified by the transmitter, Tx [Ref. 
17, ROHDE, p. 14]. After transmission from the antenna most of the power takes the 
direct path to the receive antenna. Some power, however, is reflected from stationary and 
moving objects such as buildings and aircraft. This tends to add or subtract from the 
direct path causing what is known as multipath interference noise[Ref. 18, SKOLNIK, p. 
18.46]. This invention removes the multipath noise from the received signal. 

The signal at the receive antenna is amplified and demodulated by the receiver, Rx 
[Ref. 17, ROHDE, p. 361 ][Ref. 19, VAN VALKENBURG, p. 23-6]. It is then passed 
through a presampling filter, sampled and converted from analog to digital [Ref. 19, VAN 
VALKENBURG, p. 28-7]. These real digital samples have a Fourier transform the real 
part of which is symmetric about zero frequency and the imaginary part of which is 
antisymmetric about zero frequency. The real samples are first frequency shifted by 
multiplying by e~ j2nH where t is time, F = jjj and y D is the sampling frequency [Ref. 
14, PAPOULIS, p. 15]. These complex samples are then passed through a low pass 
digital filter with a maximum frequency of F [Ref. 16, RABINER, p. 100] [Ref. 13, 
OPPENHEIM, p. 261]. The output of the low pass filter is resampled at a frequency of 
y , i.e., every other sample, to obtain the complex samples of the received signal. 



The received signal, after demodulation, analog to digital conversion and conversion to 
complex samples, can be modeled as 

N M 
h=1 m--M 

where 

s i = original signal 

/ = sample index 

D = sampling time interval 

n = time delay index 

/ = Doppler shift frequency interval 

m = Doppler shift index 

a mn = Multipath amplitude coefficient 

j = unit imaginary number 

e = base of natural logarithms 

7t = pi 

E = e j2nft 
t = iD 

r| . = white noise 
From the Figure, the residual is 

l N M 

n = q> - Z w Ai-i ~ Z Z V/-i.^ w 

where the w t are whitening filter coefficients and the x nm are the delay-Doppler filter 
coefficients. This equation corresponds to that of Box [Ref. 11, BOX, p. 498]. One can 
minimize the root mean squared (r.m.s.) residual 

over 

by nonliear least squares [Ref. 11, BOX, p. 504, p. 505] [Ref. 12, MARQUARDT, p. 
431]. 
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To show that this removes the Multipath and restores the original signal, define the 
whitened signal as 

i 

Pi = ^ ~ X w i s i~i 

and the modified noise as 

L 

Then, the residual is 

N M N M ■ 

n=\ m--M h=1 m=-M 

The root mean square residual is minimized when 
and 

Note that only delayed samples of the whitened signal p i , i.e., p._ n , n * 0, are 
subtracted from /?. , so that it is not possible to cancel the whitened signal p. since it is 
uncorrected with delayed versions of itself 

Adding the above two equations which describe the action of the whitening filter, one 
obtains 

L 

r i * Qi + )- X w i fe-/ + ) 
/=i 

Solving this equation, one gets 

L 

s i+t\i *'/ + Z w /(*m+*1/-/) 
/=i 

which, from the Figure, gives the output of the unwhitening filter. This is the desired 
result which restores the original signal without multipath. The x mn imply target 
differential range, range rate and magnitude. 



A simple test of this concept was made with the computer program shown in the 
Appendix. The Marquardt method for nonlinear least squares [Ref. 1 1, BOX, p. 504, p. 
505] [Ref. 12, MARQUARDT, p. 431] was programmed in the APL language in program 
NLS. The equation for the residual was programmed in DELRES with 
x nm = 0 when m*0 
and 
w, =0 
If 

w, * 0 

this corresponds to the equation of Box [Ref. 1 1, BOX, p. 498] and is a general recursive 
digital filter [Ref. 15, PAPOULIS, p. 45]. 

The original signal s,. was generated with a Gaussian random number generator in 
program NORM with root mean square (r.m.s.) of 1 and / = 100000 . This corresponds to 
one second of a 50000 Hz bandwidth signal represented by 100000 real samples. The 
data is white, zero mean, Gaussian with r.m.s. = 1. 

The received signal q i was generated with 

a o\ - a o3 - 0-1 9 otherwise = 0 

and 

1,=0 

This corresponds to two multipaths both of magnitude 0.1, one is delayed 1 sample, the 
other is delayed 3 samples. The white signal itself provides an effective noise which 
perturbs the estimates. 

The process converged in three iterations, the final r.m.s. residual = 1.002370568. The 
three coefficients x 0] ,x 02 ,x 03 are shown to be approximately = 0.1, 0.0, 0.1 matching 
a o\ > a n 5 ^o3 w fth r.m.s. errors of about 0.003. The mean residual was about 0.0087. The 
original signal to multipath noise ratio was 17 dB. After processing, the two multipaths 
were canceled to a signal to multipath noise ratio of 51 dB. 

The filter weight x mn corresponding to a moving target is a measure of the differential 
range nDc and range rate mfc where c is the speed of light. The system, with a common 
antenna, can be duplicated to receive two signals each from each of two widely separated 
transmitters. Geometric triangulation [Ref. 18, SKOLNIK, p.25.5, p.25.13] can be used to 
measure the two dimensional position and velocity of the target. Triplication can be used 
for three dimensional position and velocity. 

Another use of the system is to detect the angle of moving targets. If multiple antennas 
are provided, each connected with a system like that shown in the Figure, the angle of a 
moving target causing a multipath reflection to the receiving antennas can be measured 
by using the relative phases of the corresponding delay and Doppler complex coefficients 
across the several antennas [Ref. 18, SKOLNIK, p.3.34, p.3.35]. 
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The receiving antennas, receivers and processing system can be placed in a surveillance 
aircraft the position and velocity of which is obtained by an accurate navigation system 
such as GPS. Objects on the ground, interferers and targets causing multipath reflections 
can be processed by the system. Target position and velocity can be obtained by adding 
the position and velocity of the surveillance aircraft to the measured position and velocity 
of the target. 
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APPENDIX 
Computer Programs and Sample Run 



V NORM ZD IV 

[0] X+NORM N;Z;T;R 

[ 1 ] a ; ; 

[2] ft 961204.1712 

[3] A GENERATE N NORMAL RANDOM VARIABLES* MEAN 0, SIGMA 1 
[4] a DM IS SEED 
[5] ZH-2 f rjT+2 

[6] Z+Zp( * 10 000 )x? (x/Z)p 10000 

[7] 2Mo2)x2[l;] 

[8] i?+(-2x©2[2;] )*0.5 

[9] JMT+,(1 lo.xl?)x2 io.of 

V 1996-12-10 11.00.01 (GMT-*) 



VDELRES101V 

[0] R+Y DELRES X;N;I 

[I] A ;; 

[2] A 000323.2241 

[3] a DELAY RESIDUALS 

[4] N+pX 

[5] i?^oxr 

[6] RL\Nl+YliNl 

[7] 1+0 

[8] NEXTI :I>J+1 

[9] i?[J7+J]+r[tf+I]-J?[#+I- iJVD+ . xJST 

[10] +NEXTI* \I<(pY)~N 

[II] R+N+R 

V 2000-03-25 22.51.47 (GMT-*) 



VNLSID1V 

[0] JTLS RES;X;&X;R;PiJiARiA;BiDiEMiE;IiEE;AA 

[I] A ;i?£S; 

[2] A 800225.0914 800820 000321.1532 

[3] A NONLINEAR LEAST SQUARES 

[4] a Cff^i?. P£<72\ 15 RESIDUAL FUNCTION NAME 

[5] JMDCl;] 

[6] LX+XD12\1 

[ 7 ] EE+ » ' 

[8] 1+0 

[9] BBXTIiI+I+1 

[10] R+*.RES 9 • JP 

[II] A (5 it B)+±RES,'(X A*)' 

[12] A EE+EE, (*/ + /(J?*2), [0.5]1)*0.5 

[13] (*/ + /(i?*2),[0.5]l)*0.5 

[14] a ANALYTIC RESIDUAL PARTI ALS 

[15] a P+±(RES 9 'P'),' JP 

[16] n+ANPAR 

[17] a CALCULATE R 9 A AND B 

[18] a (5 4 B)+stRES t nX LX)' 

[19] a +CALCRAB 

[20] a NUMERIC RESIDUAL PARTI ALS 

[21] P+" 

[22] cT+0 

[23] NEXTJtJ+J+l 

[24] AJM±Mfl,« A-+A^[J]x e T=ipA'' 
[25] P+P,AJ?TA*[tn 
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[26] +NEXTJ*\J<pX 

[27] P^sK(pJn,(pPHp*)pP 

[28] ANPAR:A+($P)+.xP 

[29] B«-($P)+.xJ? 

[30] CALCRAB : 

[31] IH-*(((0 l + p4)pi4)[;l] + l£-10)*0.5 

[32] 44«-Ux2?o.x2?) + l£-10x( p 4)p 1,0x4 

[33] Jir*Jir-2)x(2)xfl)ffli|il 

[34] +NBXTI*xI<3 

[35] R+*.RE8, % JP 

[36] a (R A B )+&RES 9 % (X kX) % 

[37] + [0. 5 ]1 

[38] + 2 ),[0. 5] 1)*0. 5 

[39] XE+(X 9 EM) f [0. 5]£x(( (0 l + p44)p(2><> . x£ >x&44 ) [ ; 1 ] , 1 ) * 0 . 5 

[40] p EE 9 E 

[41] E 

[42] R » ' 

[43] fl AT£ 

[44] ft +NEXTI*xI<20 

[45] ft i?i?+J? 

V 2000-03-29 23.00.55 (GMT-H) 

VMS [ 04 5] 

[45] ft i?J?+i? 
[45] 

[46] 7 



pX5+N0RM 100000 

100000 

0^S^(*/ + /((Jir5)*2),[.5]l)*.5 
1.002357739 

pA"X5-«-jSr5[3+\99997]+.lxJSr5[( i99997 )o.+-l+ X 4] + .xl 0 10 

99997 

CMT<-(4/+/( (AT5-3+;5f5)*2), [. 5]1 )*.5 
0.1418755036 

20xlO©S*tf 
16.98230664 

(2 3p0 0 0 .001 .001 .OODNLS 'XX5 DELRES ' 
1.011856929 
1.002412348 
1.00237057 
1.002370568 

0.09896411424 0.001084301862 0.09780701046 0.00874503690 7 
0.003147221067 0.00 316310 5 50 5 0.0031469113 59 1.002370568 

□+J\T+( */ + /(( j?f?-6+AT5)*2), [.5]!)*. 5 
0.002841730873 

20xlO©S*ilT 
50.94879605 
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